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Laws of Large Numbers of Subgraphs in 
Directed Random Geometric Networks 

Yilun Shan^ 
Abstract 

' Given independent random points = {Xi,--- ,X„} in R^, drawn according 

CD ! to some probability density function / on M^, and a cutoff r„ > we construct 

' a random geometric digraph G(A'„, 3^„, r-n) with vertex set Xn- Each vertex Xi is 

assigned uniformly at random a sector Si, of central angle a with inclination Yi, in a 
Qh' circle of radius rn (with vertex Xi as the origin). An arc is present from Xi to Xj, if 

I Xj falls in S'i. We also introduce another random geometric digraph G(A'„,72.„) with 

vertex set Xn = {^i, • • • , A„} in R'^, d > 1 and an arc present from Xi to Xj if | |Ai — 
Xj ll < i?„^j. Here {i?„,i}i>i are i.i.d. random variables and we may take an arbitrary 
" ' ' norm 1 1 • 1 1 . In this paper we investigate two kinds of small subgraphs — induced and 

il isolated — in the above two directed networks, which contribute to understanding the 

■ local topology of many spatial networks, such as wireless communication networks. 
I We give some strong laws of large numbers of subgraph counts thus extending those 

• ■ results of Penrose [Random Geometric Graphs, Oxford University Press, 2003]. 

^ ■ 

> . 

■ Keywords. Random scaled sector graph; Random geometric graph; Law of large 
I number; Azuma's inequality 

^ ! 1. Introduction 

QQ . In the last decade there has been a resurgence of interest in the analysis of random 

[ geometric graphs (RGGs) particularly in the context of ad hoc wireless networks. An 

OA ' elegant written tutorial of random geometric graph theory is available in p^J, and the 

. paper |8j is a more recent survey emphasizing wireless networks. An RGG is usually 

O I constructed as follows. Let || ■ || be some norm on M.'^, and r„ be a real sequence. Let Xn = 

' {Xi, X2, ■ ■ ■ ,Xn}, {Xi} are i.i.d. random d-vectors in M"^ having a common probability 

^ ■ density function /. We denote by G(A'„,r„) the graph with vertex set Xn and with an 

• ^ , edge XiXj if and only if \\Xi — Xj\\ < r„ for i ^ j. Note that G{Xn-,Tn) is isotropic 

^ ' and thus undirected, which is less appropriate in many practical applications such as 

■ wireless sensor networks. The issue of small subgraph counts are dealt with here for 
two kinds of directed models of RGG (for formal definitions see below). In the case of 
wireless networks, we are sometimes interested to know the number of a desirable local 
configuration involving a small number of transmitters and receivers. This is especially 
true if there are a small number of nodes with special capabilities, e.g. data collection 
centers in sensor network process the data collected by the beacon nodes that help in self- 
organization of the network. The small subgraph counts are also of independent interest 
in random graph theory in various guises. The concerned results on small subgraph for 
classical Erdos-Renyi random graphs are discussed in detail in |2j(chap. 4) and ^(chap. 
3), and for asymptotic results in random geometric graphs, see [13] (chap. 2) and [5], while 
for exact formulae treated in the circumstance of wireless network, see (15]. 

In this short paper, we extend the method of Penrose [13] and establish some strong 
laws of large numbers of small subgraphs in directed geometric networks in some limiting 
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regimes. We now define two random geometric digraphs to be used in this work. The 
archetype (with uniformly distributed points in [0, 1]^) of the first one has been proposed 
in [B], cahed "random scaled sector graph", to model the "Small Dust" sensor networks 
using optical communication. Some graph theoretic properties have been addressed for 
this model, see e.g. [01 [71 [121 UHl Hj j mainly using combinatorial techniques. 

Definition 1. Let || • || be Euclidean norm equipped on M?. Let a € (0,2tx]. Let Xn = 
{Xi,X2,--- he i.i.d. random vectors with a common density function / : ^ M. 

Let yn = {Yi,Y2, ■ ■ ■ ,Yn} be i.i.d. random variables, uniformly distributed on [0,2ti). 
Associate every point Xi G a sector, which is centered at Xi, with radius amplitude 
a and elevation Yi with respect to the horizontal direction anticlockwise. This sector is 
denoted as S{Xi,Yi,rn). We denote by 3^„, r„) the digraph with vertex set Xn, and 

with an arc {Xi,Xj), i ^ j, presents if and only if Xj € S{Xi,Yi,rn)- 

For technical reasons we always assume r„ — > as n ^ oo. We also assume that / is 
bounded and a.s. continuous throughout the paper. We mention that the above assump- 
tions imposed on / is rather mild; in fact typical distributions such as normal distribution 
and / = l[o,i]rf are clearly satisfied. Now we introduce another model G{Xn,TZn) with 
random cutoffs motivated by Boolean model in continuum percolation [11]. For the sake 
of convenience, we still choose to employ the signs Xn, f and || • || with a little ambiguous 
(see below), however, the right meaning will be clear in the context and no confusion will 
be incurred. 

Definition 2. Let || • || be any norm equipped on W^, d> 1. Let Xn = {Xi,X2, ■ • • ,Xn} 
be i.i.d. random vectors with a common density function / : M'^ ^ M. Let Rn be a positive 
random variable with probability distribution function Qn and density function qn- For each 
point Xi, we associate a ball B{Xi, Rn^i) with radius Rn,i, centered at Xi, independent of 
other points. {Rn,i}^^i are independent copies of Rn and setTZn ■= {Rn,i, Rn,2, ■ ■ • ,Rn,n}- 
We denote by G{Xn,'JZn) the digraph with vertex set Xn, and with an arc {Xi,Xj) origi- 
nating from Xi and terminating in Xj if and only if Xj € B{Xi, Rn,i)- 

As usual, we shall impose a certain decaying condition on Here we assume ERn = 

r'^dQn{r) — > as n oo, throughout the paper. We will investigate two kinds of 
subgraphs in the above two models; one is induced subgraph and the other is isolated 
subgraph. Induced subgraph is defined in its usual meaning, see e.g.[T3]. Suppose G is a 
digraph, a subgraph H is isolated in G if is an induced subgraph and there are no arcs 
leave H. If G is undirected, then a connected isolated subgraph of G is just a component. 

Before going further we will need some other definitions. Given a finite set X in M*^, 
let card(Af) denote the number of points in X and let [ • [ be the d-dimensional Lebesgue 
measure, which is easy to discriminate from the sign for absolute value in the context. In 
the rest of the paper, let /max := sup{t : \{f{x) > t}\ > 0} be the essential supremum of 
the probability density function /. As mentioned before, we assume /max < oo. For a set 
A € M'^, let X{A) denote the number of points of X located in A. Denote D{0, 1) as the 
unit disk in M^, then \D{0, 1)| = ti and also set 6 := \B{0, 1)| w.r.t some given norm. Let 
C, G' etc. be various positive constants, and the values may change from line to line. 

The rest of this paper is organized as follows. Section 2 contains the statement of main 
results and proofs are provided in Section 3. We finally draw conclusions in Section 4. 

2. Statement of main results 

For A; G N, let T be a fixed connected graph on k vertices. We say T is feasible if 
either P{G{Xk,yk,r) ^ T) > for some r > or P(G(A'fc,7^„) ^ T) > for some 
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{rn,i,i"n,2, ■ ■ ■ , i"n,k}- Let Hn,T and Hn^T be the number of induced subgraphs and isolated 
subgraphs of 3^„, r„) isomorphic to T (T-subgraphs) respectively. Likewise, let Gn,T 

and Gn,T be the number of induced and isolated T-subgraphs of G{Xn-,'R-n) respectively. 
For a finite set A' C 1^2 ^nd a point y € [0, 2nY''^''^^^\ we define indicator random variables 
hT{X,y) ■■= \G{xy,i)?^T\ and hn,T{X,y) := l[G{x,y,r„)^T]- For a finite set ^ C R'^ and 

TZn = {Rn,l, Rn,2,- ' ' , Rn,ca.id{X)} , We define gn,T{X.,TZn) ■= '^[G{X,'R.n)^T]- 

The basic tool we shall need in the proofs is the following Azuma's inequality; we refer 
the reader to [1] for a proof and a wealth of materials regarding that topic. 

Lemma 1 . Suppose Mi , M2 , Mn is a martingale with corresponding martingale dif- 
ference sequence Di, D2, • ■ ■ ,Dn, where Di := Mj — Mj_i [2 < i < n) and Di := 
Ml — EMi . Then for any a > 0, we have 
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P(| a| > a) < 2exp ( - n° , 



where := inf{& : -P(|A| < &) = !}• 

In the sequel we sometimes use a generalized version with "tolerance" of Azuma's 
inequality [3]: 

Lemma Suppose Mi, M2, • • • , M„ is a martingale with corresponding martingale 

difference sequence Di,D2, — ' ,Dn, where Di := Mi — Mj_i (2 < i < n) and Di := 
Ml - E'Mi . Then for any a,b > 0, 

p(|Ea|>«).2.p(-^).(i.?i«)Ep(iai>.). 

i=l i=l 

By the notations defined in the beginning of this section, we have the following strong 
laws of large numbers under various regimes: 

Theorem 1. Suppose T is a connected feasible graph of order k, k £ N. Let nr'^ ^ A G 
(0,00). Then 

lim n~^HnT = k~^ / ipT{Xf{x))f{x)dx a.s. 

where 

/'2ti /»2ti r r 



(pT{t) :-- 



k k-1 

.g-t|S(o,x2,- ,xk,yi,- ,yk)\^x2 ■ • • dxfcdyi ■■■dyk k>2; 
6 2 A; = 1, 



and S{xi, ■■■ ,Xfc,yi, • • • ,yk) ■= ^i=iS{xi,yi, 1). 

Theorem 2. Let T be a single vertex. Suppose nR^ — > A S [0, 00) in probability. Then 



lim n-^GnT=f e-^^^(^)/(^)dx a.s. 



Theorem 3. For k > 2, let T be a connected feasible graph of order k. Suppose nr"^ 
andlnn = o{n?'^~^r'^^ as 00. Then 

lim n''''r:^^^'^^^Wn,T = a.s. 
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where 

t^T ■■= TTT^-TT I •••/ / ••• / hTiO,Z2r-- ,Zk,yi,--- ,yk) 




k 

■ /^(x)dxdz2 • • • dzfcdyi • • • dyk- 

Theorem 4. For k > 2, let T be a connected feasible graph of order k. Suppose f has 
bounded support set (denoted by suppf). Suppose r„ ^ and there is a constant 6 > 
such that liminf n^^~-'^~'^r.^^^ > 0. Then 

\mi n~^r~'^^^~''^HnT = IJ-T a.s. 

n— >oo ' 

3. Proofs 

To prove our main theorems, we first derive several asymptotic resuhs (Propositions 
1~4) for the means of subgraph counts -ffn,T, -^n,T and Gn,T- 

Proposition 1. For k > 2, let T be a connected feasible graph of order k. Suppose 
r„ — > as n ^ oo. Then 

hm n-^r-'^^^-^^EHnT = l^T 

n—foo ' 

where ht is defined in Section 2. 

Proof. From the linearity of expectation, EHn^T = (fc)-E'^n,r('^fc) 3^A:)- Thereby 



2n f2n 



•/''(xi)dxi • • -dxfcdyi ■■■dyk 

+7^T~Tr(?) / •■■ / / ■■■ / K,t{xi,--- ,Xk,yi,--- ,yk) 



(2ti)'= 



Jo JR2 



■ ( n ~ f^i^i) j dxi • • • dxkdyi ■■■dyk 



Let xi = X and Xi = xi + rnZi for 2 < i < /c, then the first term on the right hand side of 
([1]) equals 



TTT-TT ( ? Vn ■■■ / / ■■■ / K,T{x,x + rnZ2r-- ,x + rnZk,yir-- ,yk) 

■ f''{x)dxdz2 ■ ■ ■ dzkdyi ■ ■ ■ dyk- 

We have hn^rix, x+rnZ2, ■ ■ ■ ,x+r„Zfc,yi, • • • ,yk) = /it(0, ^2, • • • • • • ,yk)- Since / is 

bounded and T is a connected graph, /^(x) and hxiO, Z2,^ ■ ■ , Zk,yi, • • • , yk) are integrable 
on M? and (M^)'^^^ x [0,2ti)'^ respectively. Hence, the first term on the right hand side of 
([1]) tends to n^r"^^ by the dominated convergence theorem. 

Set rin{xi) := • • • r^^^^^^"^) | /(^i) - /^"Ha^i)|dx2 • ■ ■ dx^. It's 

easy to see that the absolute value of the second term on the right hand side of ([I]) 
multiplied by n~^rn'^'^^ is bounded by j^2f]n{xi)f{xi)dxi. If xi is a continuous point 
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of /, then Tjnixi) — > as n — ^ cx) by the mean value theorem for integrals. Therefore, by 
the dominated convergence theorem and the assumption of almost everywhere continuity 
of /, we have J^2 Vn{xi)f{xi)dxi — > 0. This concludes the proof. □ 

Proposition 2. For k > 2, let T he a connected feasible graph of order k. Suppose 
nr'^ — > as n ^ oo. Then 

lim n-V-2(*=-i)E^,,T = //T 

n-+oo 

where fj,T is defined in Section 2. 

Proof. Let hn,T{X,y) := 'i-[G{x,y,rn) is isolated and ^T]- Therefore, 

EHn,T = ij^]Ehn,T{^k,yk) 



(:)■ 



\Ehn,TiXk,yk) ■ P{G{Xk,yk,rn) is isolated|G(A'fc,3'fe,r„) ^ T) 
:= EHn,T-Pi 

Since T is a connected graph and by the assumed asymptotic behavior of r„, we get as 
n — oo, 

1 > Pi > (1 - PiXi e L»(0, krnW-" > (1 - f.n.Mkrnfr'" ^ 1- 

Consequently, EHn,T = (1 + o{l))EHn^T- By using Proposition 1, we complete the proof. 
□ 

Proposition 3. For k € N, let T be a connected feasible graph of order k. Suppose 
nr"^ — >■ A G (0, oo) as n ^ oo. Then 



lim n-^EHn,T = k-^ [ ^T(A/(x))/(x)da 



where 'i^ri ) is defined in Section 2. 

Proof. By the definition of hn,T and similar with the beginning of the proof of Proposition 
1, we have 

1 fn^ '•^^ r2n 



n{2n)>' \kj Jo Jo Jm.^ Jr^ 



• (l - / fix)dx) f''{xi)dxi ■ ■ ■ dxkdyi ■■■dyk 

— T^S^frl / ■■■ / I "I hn,Tixi,--- ,Xk,yi,--- ,yk) 

n{2nf \kj Jo Jo 7m2 7m2 

■(i- [ f{x)dxy~' 

k . 

n /(^i) - ) da^i • ■ ■ d^fedyi ■■■dyk (2) 

i=i ^ 
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Let Xi = xi + rnZi for 2 < i < k, then the first term on the right hand side of ([2]) tends to 



1 / 



/o \k\^]^n'' ••• / / / K,Tixi,Xi+rnZ2,--- ,Xi+rnZk,yi,--- ,yk) 

n{2n)>' \kj Jq Jq J^2 J^-z 

■ e^""'')^'' (i-4(.i,.i,.„)uut2S(-i+^n....„.n) fk^xi)dxidz2 ■ ■ ■ dzfcdyi • • • dyfc 



which is further asymptotic to 



/c!(2tx)'^ 







hT{0,Z2,--- ,zk,yir-- ,yk) 

2 



by using the mean value theorem for integrals and the dominated convergence theorem. 
Thereby, the first term on the right hand side of ([2]) tends to y^T{^f{xi))f{xi)dxi. 

On the other hand, note that n~^{2^ < Cvn'^^^ for some positive constant C, 
then the absolute value of the second term on the right hand side of ([2|) is bounded by 
C/jj2 /(xi)7?„(2;i)dxi, where r/„(xi) is defined by r/„(xi) = Jo^ccukm) ' ' ' /-D(xi,fcr„) ^n'^^'''^^ 
■ 1 11^=2 /(-^i) ~ f^~^(xi)\dx2---dxk- By the mean value theorem for integrals, rin{xi) 
tends to if xi is a continuous point of /, then by the dominated convergence theorem, 
as in the proof of Proposition 1, J^2 /(xi)r/„(xi)dxi ^ as n — > oo. Thus the proof is 
completed. □ 

Proposition 4. For k = 1, let T be a graph of order k, that is, T is a single point. 
Suppose nR^ ^ A S [0, oo) in probability, as n ^ cx). Then 

lim n-^EGnT=f e'f^^^'^^f{x)dx 
Proof. By the definition of gn,T, 
n~^EGn,T = I I an,T{xi,rn,i)(l- I /(x)dx) /(xi)g„(r„,i)dxidr„,i 

Jo JR"^ ^ JB(xi,rni) ' 



oo 



JR'^ ^ JB{xi,r,^^i) 



' B(xi,rn.i) 
\ n— 1 

/(x)dx /(xi)g„(r„,i)dxidr„,i 



e-"/(-i)e<,i/(xi)g„K,i)dxidr„,i. 

Hence using the dominated convergence theorem and the assumption of Rn, the above 
expression tends to J^a e--^(^i)^^/(xi)dxi as n tends to infinity, which concludes the proof. 
□ 

Now we are in position to prove our strong laws of large numbers. 

Proof of Theorem 1. In order to use Lemma 1, we shall first define a filtration. Let 
To = {0, ^i} be the trivial u-field, and J^i = a{{Xi, Yi), • • • , {Xi,Yi)} for i > 1. Define a 
martingale difference sequence as Dn,i ■= E{Hn,T\J^i) — E{Hn,T\J^i-i), therefore we can 
write Hn,T — EHn,T = Y17=i ^n,i- Let j, denote the number of isolated T-subgraphs in 
G{Xn+i\{Xi},yn+i\{Yi},rn) and H^^irp denote the number of isolated T-subgraphs in 
G{Xn+i,yn+i,f'n). Thus we have Dn^i = E{Hn,T — H'!^rp\Ti). Adding a point to a finite 
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set in can cause the number of isolated T-subgraphs to increase by at most 1, and 
can cause it to decrease by less than a constant M (here d = 2, so we may take M = 6), 
thereby we get 

\Hn,T - H'^^tI < \Hn,T - H!^+1,t\ + \H'n+l,T - Hn,T\ 

< (M + 1) + (M + 1) = 2(M + 1) 

and then Dn^i < 2(M + 1). Now for any e > 0, by Lemma 1, we have P{\Hn,T — EHn,T\ > 
en) < 2exp ( — g^^^pjp); which is summable in n. The result follows by Borel-Cantelli 
Lemma and Proposition 3. □ 

Proof of Theorem 2. The proof parallels to that of Theorem 1. Define a filtra- 
tion: J^nfl = {0,^^} and J^n,i = cr{{Xi, Rn,i), ■ ■ ■ , {Xi, Rn^i)} for i > 1. A martingale 
difference sequence is defined by Dn,i = E{Gn,T\^n,i) — E{Gn,T\^n,i~i), and then we 
have Gn,T — EGn,T = '127=1 ^n,i- Let G^x denote the number of isolated vertices in 
G(A'„+l\{Xi},7^„ U Rn,n+i\{Rn,i})- ^^'s easy to see that Dn,i = E{Gn,T - G\^r^\Tn,i)- 
Reason similarly as in the proof of Theorem 1, there exists a constant M > depends 
only on d such that \Gn,T — G^tI — hence \Dn,i\ < M a.s. For e > 0, by Lemma 

1, P{\Gn,T — EGn,T\ > £n) < 2 exp ( — j^), which is summable in n. The result then 
follows by Borel-Cantelli Lemma and Proposition 4. □ 

Proof of Theorem 3. Proceeding along the same line as the proof of Theorem 1, by 
Lemma 1, we get 

(2 2fc— 1 4(fc— 1) \ 

which is summable in n by the assumption. The result follows by Borel-Cantelli Lemma 
and Proposition 2. □ 

In the next proof, we substitute Lemma 2 for Lemma 1 since this time the basic Azuma' 
inequality is no longer valid. 

Proof of Theorem 4. Set r] := ^^^'^-^^ and partition M? into squares {A^^i, i G N), each 
of side r„. Let 

E := {X C M^|card(<%') = n, Af(yl„_i) < n^{nrl V 1) for every An,i s.t. A^a n supp/ / 0}. 

Since Xni^An^i) ~ Bin{n, f{x)dx) is a binomial random variable, by a Chernoff bound 

(see e.g.[l3] p. 16), we have P{Xn{An,i) > n^{nr'^ V 1)) < e^"'' when n is large enough. 
Since supp/ is bounded, we have 

P{Xn ^E) < Cn^^e-''\ (3) 

where C is some positive constant. Now let's define a filtration. Let J-q = {0,1^} be the 
trivial a- field, and J^i = (t{(Xi,1i), • • • , {Xi,Y,i)} for i > 1. We have Hn,T — EHn,T = 
Y17=i with the martingale differences Dn,i ■= E{Hn,T\J^i) — E{Hn,T\J^i-i)- Let H^rp 
denote the number of induced T-subgraphs in G{Xn+i\{Xi} , yn^i\{Yi} , rn) , then Dn^i = 
E{Hn,T — Hl^rp\J^i). Set an event En,i = {Xn G E, Xn+i\{Xi} G E}, hence we may derive 
\Hn,T - H'!^^t\ • ^E„,i < C ■ [n'^{nr'^ V 1)]^"^ for some constant C, and \Hn,T - -f^^,rl ^ 
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since changing the position of one point in a configuration can only at most affect the 
subgraphs constructed by points in surrounding nine squares. Consequently, 

\Dn,i\ < E{\Hn,T - HI^,t\ • '^E^JJ^i) + E{\Hn,T " -f^n,rl ' l^^^J-^i) 

< C-[vP{nrlyl)f-'+n^P{El^,\T,) (4) 
Define an event Fn^i := jjJFj) < n^*^}, thus 

P(F;,) < n'E{P{El^,\T.)) = n'P{E^n,) < Cn'+'-^e-^', 

by Markov's inequality and ([3]). Note by ([5]) when occurs, |-Dn,i| ^ C ■ [n"^ {nr^W . 
Therefore, by using Lemma 2, we obtain for any e > 

/ 2 2k 4(fc— 1) \ 

which is summable in n. The result then follows by Borel-Cantelli Lemma and Proposition 
1. □ 

4. Conclusions and future work 

We have derived various strong laws of large numbers of subgraph counts in two types 
of directed random geometric networks. These models of directed graphs are natural 
generalizations of standard geometric graphs applicable to analysis of a variety of spatial 
networks exemplified as wireless communication networks and sensor networks. Our results 
obtained for the second model is limited and a very natural question would be to ask what 
happens for other sorts of small subgraphs. 
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